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Abstract 

Events with a final-state proton carrying a large fraction of the proton beam mo- 
mentum, xl > 0.6, and the square of the transverse momentum p\ < 0.5 GeV 2 , 
have been studied in e + p collisions with the ZEUS detector at HERA. Data with 
different photon virtualities were used: Q 2 < 0.02 GeV 2 , 0.1 < Q 2 < 0.74 GeV 2 
and 3 < Q 2 < 254 GeV 2 , corresponding to integrated luminosities of 0.9, 1.85 
and 3.38 pb -1 . The cross sections are given as a function of xl, Pt, Q 2 and 
the Bjorken scaling variable, x. The ratio of the cross section for leading proton 
production to the inclusive e + p cross section shows only a mild dependence on 
Q 2 and on x. In the region 0.6 < Xl < 0.97, the relative yield of protons is only 
a weak function of xl. 
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1 Introduction 



Events with a final-state proton carrying a large fraction of the available energy but a 
small transverse momentum have been studied in detail in high-energy hadron-proton 
collisions [1]. The cross section for such leading proton events shows a peak for values of 
the final-state proton momentum close to the maximum kinematically allowed value, the 
so-called diffractive peak. For lower momenta, the cross section is lower and the fraction 
of events with a leading proton is approximately independent of the energy and type of 
the incoming particle. This characteristic behaviour has led to studies of the associated 
event in terms of the effective energy available for hadronisation [2,3]. More recently, 
events with neutrons or protons carrying a large fraction of the proton-beam momentum 
have also been measured in positron-proton (e + p) scattering at HERA [4-8]. 

The study of semi-inclusive rates in hadron-nucleon collisions indicates that particle pro- 
duction from the target nucleon is independent of the type of the incident hadron, a prop- 
erty known as vertex factorisation [1]. This has been studied for instance by comparing 
semi-inclusive rates, normalised to the respective total cross sections, for the production 
of particles in the fragmentation region of the target nucleon. The hadronic data [9] also 
show that, in the high-energy limit, the momentum distribution of the particles from the 
fragmentation of the target hadron is independent of the energy of the incoming parti- 
cle. These characteristics have not yet been extensively studied for baryon production in 
electron-proton collisions. 

Electroproduction of leading baryons is also interesting in other respects. The virtual 
photon mediating the interaction, in the reference frame in which the proton is at rest, 
fluctuates into a vector-meson-like object (the vector dominance model, VDM [10]). The 
transverse size of this projectile can be varied by changing the virtuality, Q 2 , of the 
photon. Real photons (Q 2 = 0) have hadronic size, while, as Q 2 increases, the photon 
size decreases. It is thus possible to experiment with a projectile of varying size. This 
may lead, for instance, to different absorptive rescattering of the produced baryon as Q 2 
changes [11], and hence to a breaking of vertex factorisation. 

Several mechanisms have been suggested to explain the hadroproduction or electropro- 
duction of leading protons. None of them are as yet amenable to calculations based on 
perturbative quantum chromodynamics (pQCD). This is, in part, a consequence of the 
fact that the small px of the leading proton necessitates a non-perturbative approach. 
Some models [12-16] are based on the Regge formalism, with leading proton production 
occurring through t-channel exchanges, both isoscalar and isovector, notably Pomerons, 
Reggeons and pions. These exchanges mediate the interaction between the proton and 
the fluctuations of the virtual photon. The relative contribution of the different exchanges 
varies as a function of the momentum and type of the scattered baryon: for leading pro- 
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tons, Pomeron exchange dominates in the diffractive-peak region with Reggeon and pion 
exchanges contributing for lower outgoing-proton momenta. Other theoretical models re- 
tain quarks and gluons as fundamental entities, but add non-perturbative elements, such 
as soft-colour interactions [17]. The concept of fracture functions also offers a general 
theoretical framework for a QCD-based study of leading baryon physics [18]. 

This paper reports studies of leading proton production in e + p collisions, e + p — > e + Xp, 
emphasizing the non-diffractive region. This complements the recent ZEUS study of 
leading neutrons [8]. High-energy protons with low transverse momentum carrying at 
least 60% of the incoming-proton momentum were measured in the ZEUS leading proton 
spectrometer (LPS) [19]. The fraction of such events with a large rapidity gap in the 
forward region is presented. The longitudinal- and transverse-momentum spectra are 
studied for different photon virtualities, from quasi-real photoproduction (Q 2 ^0.02 GeV 2 ) 
to Q 2 = 254 GeV 2 . The dependence of the cross section for the production of leading 
protons on the deep inelastic scattering (DIS) variables x and Q 2 is also measured and 
compared to that for the inclusive reaction e + p — > e + X . The results are discussed in the 
context of vertex factorisation and particle-exchange models. Finally, the properties of 
events with a leading proton and two jets are presented. 

2 Experimental set-up 

The measurements were performed at the DESY ep collider HERA using the ZEUS detec- 
tor. In 1994 and 1995, HERA operated at a proton energy E p = 820 GeV and a positron 
energy E e = 27.5 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [20]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked by the central tracking detector (CTD), which operates in a 
magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists of 72 
cylindrical drift-chamber layers, organised in nine superlayers covering the polar-angle 1 
region 15° < 6 < 164°. The relative transverse-momentum resolution for full-length tracks 
is a{pt)/p t = 0.0058pt © 0.0065 © 0.0014/p t , with p t in GeV [21]. 

The high-resolution uranium-scintillator calorimeter (CAL) [22] consists of three parts: 
the forward (FCAL), barrel (BCAL) and rear (RCAL) calorimeters. Each part is subdi- 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — ln(tan §), where the polar angle, 8, is measured with respect to the proton beam 
direction. 
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vided transversely into towers and longitudinally into one electromagnetic section (EMC) 
and either one (in RCAL) or two (in FCAL and BCAL) hadronic sections (HAC). The rela- 
tive CAL energy resolutions are c(E) / E = 0.18/ \[E for electrons and cr(E) / E = 0.35/ \[E 
for hadrons (E in GeV). 

A lead-scintillator calorimeter (LUMI-e) at Z = —35 m [23], with a relative energy reso- 
lution oia(E)/E = 0.18 j\fE (E in GeV), was used to tag events with positrons scattered 
through angles up to about 5 mrad, and to measure the scattered-positron energy, E' e , over 
the range 7 < E' e < 21 GeV. These events have Q 2 < 0.02 GeV 2 and are hereafter referred 
to as the "photoproduction" sample. This sample was collected in 1994 and corresponds 
to an integrated luminosity of 0.898 ± 0.014 pb" 1 . 

Low-Q 2 events (0.1 < Q 2 < 0.74 GeV 2 ) were tagged by requiring the identification of the 
scattered positrons in the beam pipe calorimeter (BPC) [24-27], a tungsten-scintillator 
sampling calorimeter, located close to the beam pipe, 3 m downstream of the interaction 
point in the positron beam direction. This low-Q 2 sample, hereafter referred to as the 
"BPC sample", has an integrated luminosity of 1.85 ±0.02 pb _1 and was collected in 1995. 

For higher-Q 2 events (Q 2 > 3 GeV 2 ), the impact position on the CAL surface of the 
scattered positron was determined with the small-angle rear tracking detector (SRTD) [28] 
or the CAL. The SRTD is attached to the front face of the RCAL and consists of two planes 
of scintillator strips, 1 cm wide and 0.5 cm thick, arranged in orthogonal orientations. 
Events with Q 2 > 3 GeV 2 are referred to as the "DIS sample" in the following. The 
integrated luminosity of this sample, which was collected in 1995, is 3.38 ± 0.03 pb -1 . 

The leading proton spectrometer (LPS) [19] detected charged particles scattered at small 
angles and carrying a substantial fraction of the incoming-proton momentum; these par- 
ticles remain in the beam pipe and their trajectories were measured by a system of silicon 
micro-strip detectors inserted very close (typically a few mm) to the proton beam. The 
detectors were grouped in six stations, SI to S6, placed along the beam-line in the di- 
rection of the proton beam, between Z = 20 m and Z = 90 m. The track deflections 
induced by the magnets in the proton beam-line allow a momentum analysis of the scat- 
tered proton. During data taking, the detector planes were inserted close to the beam 
by means of re-entrant pots and were retracted during beam dump and fill operations 
of the HERA machine. For the present measurements, only the stations S4, S5 and S6 
were used. The intrinsic resolution is better than 1% on the longitudinal momentum and 
5 MeV on the transverse momentum. The effective transverse-momentum resolution is, 
however, dominated by the intrinsic transverse-momentum spread of the proton beam at 
the interaction point, which was ~ 40 MeV in the horizontal plane and ~ 100 MeV in the 
vertical plane. 
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3 Kinematics and cross sections 



Figure 1 illustrates semi-inclusive leading proton production in ep collisions. Four kine- 
matic variables are needed to describe the interaction e + p — > e + Xp. They are defined in 
terms of the four-momenta of the incoming and outgoing positron, K and K', and of the 
incoming and outgoing proton, P and P', respectively. 

Two of the kinematic variables were chosen from among the Lorentz invariants used in 
inclusive studies, of which only two are independent: Q 2 = —q 2 = —{K — K') 2 , the 
virtuality of the exchanged photon; x = Q 2 / (2P ■ q) and y = q- P '/ \K • P) ~ Q 2 / (sx); and 
W 2 = (P + K — K') 2 = m 2 + Q 2 (l —x)/x, the square of the photon-proton centre-of-mass 
energy. In these equations, m p is the mass of the proton and y/s = 300 GeV is the e + p 
centre-of-mass energy. 

Two additional variables are required to describe the leading proton. They can be chosen 
as the momentum fraction carried by the outgoing proton 

P'-K 
XL= P~K 

and its transverse momentum with respect to the direction of the incoming proton, px- 
In terms of these variables, the square of the four-momentum transfer from the target 
proton is given by 

X L X L 

where the second term is the minimum kinematically allowed value of \t\ for a given xl- 
In a particle-exchange model, t is the square of the four-momentum of the exchanged 
particle. The p\ range covered by the present data, < p% < 0.5 GeV 2 , thus translates 
into < \t\ < 0.5 GeV 2 at x L = 1 and 0.2 ^ \t\ ^ 1 GeV 2 at x L = 0.6. 

The differential cross section for inclusive e + p — > e + X scattering is written in terms of 
the proton structure function, F 2 (x, Q 2 ), as 



V + — F 2 (x,Q )(! + A), 



dxdQ 2 xQ 4 V 2 

where A is a correction that takes account of photon radiation, Z° exchange, and the 
longitudinal structure function, Fl- In analogy with this, the differential cross section for 
semi-inclusive leading proton production, e + p — > e + Xp, is written as 

■ W f i-y + y 4) F 2 H4> (^ Q 2 > ^ P 2 T )(i + a lp ), (2) 



dxdQ 2 dxLdp\ xQ 4 
where A LP is the analogue of A. 
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3.1 Reconstruction of the kinematic variables 



Three samples of data were used: 

• the photoproduction sample, with the scattered positron tagged in the LUMI-e calorime- 
ter; 

• the BPC sample, with the scattered positron measured in the BPC; 

• the DIS sample, with the scattered positron detected in the CAL. 

Different methods were used for the reconstruction of the kinematic variables, Q 2 and W, 
for the three data sets. Tagging of the scattered positron in the LUMI-e calorimeter for 
photoproduction events does not allow the measurement of Q 2 event by event; however, 
the angular acceptance of the LUMI-e calorimeter limits the Q 2 range to the region 
Q 2 < 0.02 GeV 2 . For these events, W was measured from W 2 = ys, with y = (E e —E' e )/E e , 
where E' e denotes the energy of the outgoing positron. For the BPC sample, E e and the 
positron scattering angle, i? e! as measured in the BPC, were used ("electron method" [24]) 
to determine Q 2 , W, x and y. For the DIS sample, these variables were reconstructed 
using the double angle method [29]. 

For the reconstruction of the hadronic final state, X, the energy deposits in the CAL 
and the track momenta measured in the CTD for the charged particles were clustered 
into energy-flow objects (EFOs) which are assumed to correspond to particles and are 
assigned the pion mass [30,31]. The EFOs were used to reconstruct the mass, M x , of the 
hadronic final state contained in the central detector. Using the EFOs, the y variable was 
also reconstructed with the Jacquet-Blondel method [32], which uses information from 
the hadronic final state to reconstruct the event kinematics, and was denoted by y,jB- 
Furthermore, the variable 

6 = '%2(E i -p z J + E' e (l-coBd e ) 

i 

was evaluated, where £\ denotes a sum over all EFOs, excluding those assigned to the 
scattered positron, and Ei and pz,i are the energy and the longitudinal momentum of 
each EFO, respectively. For perfect resolution and fully contained events, energy and 
momentum conservation constrain 5 to be twice the positron beam energy. The angle of 
the hadronic final state (as measured in the ZEUS central detector) with respect to the 
incoming-proton direction was evaluated from 

(Eifir,i) 2 -(Ei(^-p,,t)) 2 

where the sums run over all EFOs excluding those assigned to the scattered positron. 
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The modulus of the momentum of the scattered proton, p', was measured in the LPS, 
along with its component perpendicular to the mean proton beam direction, p?. The 
variable xl was evaluated as = p'/E p . The mean direction of the incoming proton 
beam was determined for each proton fill of HERA using the reaction ep — > ep°p at 
Q 2 « [19]. 

In the following, the term "leading proton" is used to indicate a positively charged particle 
detected in the LPS. In the present measurement, Xl is restricted to values larger than 0.6. 
Charged-particle production measured at the ISR [2, 33] shows that the pion-to-proton 
ratio at xl = 0.6 is about 10% and falls rapidly for increasing values of xl- 

4 Event selection 

Photoproduction events were selected at the trigger level by requiring a coincidence be- 
tween an energy deposit of at least 5 GeV in the LUMI-e and of at least 464 MeV in 
the RCAL (excluding the towers immediately adjacent to the beam-pipe) or at least 
1250 MeV (including those towers). This requirement helps to suppress the background 
from bremsstrahlung events {ep — > e^p) characterised by having a scattered positron in 
the LUMI-e and no activity in the rest of the detector. The BPC and DIS events were 
triggered by requiring the presence of a scattered positron in the BPC and the CAL, 
respectively. No requirement was imposed on the final-state proton at the trigger level. 

Events were selected offline in three steps: first inclusive events with the scattered positron 
in the LUMI-e, the BPC or the CAL were selected; then, a track in the LPS was required; 
finally a search for jets was carried out in the hadronic final state, X. The details of the 
selection procedure are presented in the rest of this section. The following requirements 
were imposed for all samples: 

• the Z coordinate of the reconstructed vertex, if measured, in the range —50 < Z < 
50 cm; 

• the timing of the interaction, as measured by the CAL, consistent with the timing of 
an e + p bunch crossing. 

For the BPC and DIS samples, the requirement 35 < 5 < 65 GeV was also imposed 
in order to reduce the photoproduction background and to minimise the effect of the 
radiative corrections. 

The photoproduction sample [34] was selected by requiring that a positron be measured 
in the LUMI-e with energy in the range 12 < E' e < 18 GeV, corresponding to 176 < 
W < 225 GeV. Overlay events, in which some activity in the RCAL accidentally overlaps 



6 



with the scattered positron of a bremsstrahlung event (ep — > ejp) in the LUMI-e, were 
subtracted as discussed in an earlier study [35]. The subtraction was less than 3%. 

The BPC sample [24, 36] was selected by requiring a scattered positron measured in the 
BPC with E' e > 7 GeV and a photon virtuality in the range 0.1 < Q 2 < 0.74 GeV 2 . In 
addition, the requirement 0.08 < y < 0.74 was imposed, which corresponds to 85 < W < 
258 GeV. Finally, yjB > 0.05 was required, thus ensuring hadronic activity away from the 
forward direction and reducing the migration of events from low y, where the resolution 
of the electron method is poor. 

In the DIS analysis [34,36], a scattered positron with energy E' e > 10 GeV was required 
in the CAL; the photon virtuality was restricted to the interval 3 < Q 2 < 254 GeV 2 and 
W to the region 45 < W < 225 GeV. Finally, the condition yjB > 0.03 was imposed. 

The total number of events thus selected was approximately 94 000 for the photoproduc- 
tion sample, 50 000 for the BPC sample and 386 000 for the DIS sample. 

Next, high-momentum protons in the LPS were selected by requiring: 

• one track in the LPS with p 2 T < 0.5 GeV 2 and 0.6 < x L < 1.02. For x L > 0.97, a lower 
bound on p\ was also imposed: p\ > 0.073 GeV 2 . For the 1994 data, the p\ range 
was restricted to p\ < 0.04 GeV 2 . The p\ cuts and the lower limit on xl restrict the 
data to a region of well understood acceptance; 

• no candidate track was accepted if, at any point, the minimum distance of approach 
to the beam pipe, Ap ipe , was less than 0.4 mm (0.5 mm for the 1994 data). This 
cut reduced the sensitivity of the acceptance to the uncertainty on the location of the 
beam-pipe apertures; 

• events in which the reconstructed proton track passed closer than a distance A plane = 
0.2 mm to the edge of any LPS detector were rejected. This ensured that the track 
was well within the active region of the detectors; 

• the total E + p z of the event was required to be smaller than 1655 GeV. For fully 
contained events, this quantity should be equal to 2E P = 1640 GeV. This cut rejects 
random overlays of a beam-halo proton and an event satisfying the trigger and selection 
cuts applied to the non-LPS variables; 

• for xl > 0.97, M x > 2 GeV was required, where M x is the reconstructed hadronic 
mass in the CAL. This rejects contributions from exclusive production of low-mass 
vector mesons, which have Q 2 and t dependences different from those of the inclusive 
events [37]. 

After this selection, the total number of events with a good LPS track was 1834 for the 
photoproduction sample, 1697 for the BPC sample and 13335 for the DIS sample. 
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Finally, a search was performed for jets in the hadronic final state [38]. Because of the 
limited statistics of the photoproduction and BPC samples, the search was limited to 
the DIS data. The jets were reconstructed using the algorithm [39], requiring a jet 
transverse energy Et > 4 GeV in the 7*p centre-of-mass system and a jet pseudorapidity, 
in the laboratory frame, in the range —2 < r]i ct < 2.2. A sample of 225 events with exactly 
two jets was selected. 

5 Monte Carlo simulation 

Several Monte Carlo (MC) generators were used to determine the acceptance of the ap- 
paratus for events with a leading proton. The EPSOFT2.0 Monte Carlo [40-42] was used 
for the BPC data. This generator simulates diffractive processes with dissociation of the 
virtual photon, as well as non-diffractive processes. Vertex factorisation is assumed. The 
HERACLES4.6 Monte Carlo [43], which simulates initial- and final-state QED radiation, 
is interfaced to EPSOFT. Samples of DIS events were simulated with RAPGAP [44, 45] 
version 2.06/06, which incorporates meson and Pomeron exchange; it also assumes vertex 
factorisation. QED radiation was also simulated using HERACLES. Weights were as- 
signed to the events generated with EPSOFT and RAPGAP such that the reconstructed 
proton xl and p\ spectra agreed with the data. For the photoproduction data, and for 
systematic studies, events that only contain a proton with distributions in and p\ 
tuned to those of the data were generated. This simulation produced the same results for 
the LPS acceptance as EPSOFT and RAPGAP. 

All generated events were passed through the trigger-simulation package and the stan- 
dard ZEUS detector simulation, based on the GEANT 3.13 program [46]. The simulation 
includes the geometry of the beam-pipe apertures, the HERA magnets and their mag- 
netic fields. The spread of the interaction- vertex position was also simulated, as were the 
proton-beam angle with respect to the nominal direction and its dispersion at the interac- 
tion point. The simulated events were then passed through the same reconstruction and 
analysis programs as the data. 

Figures 2a)-c) compare the distributions of the reconstructed variables E' e , $ e , and tjjb 
in events generated with EPSOFT with those for the BPC data. The agreement between 
the data and the simulated distributions is good. A similarly good description of the DIS 
data by RAPGAP is shown in Figs. 2d)-f) for the variables E' e , $ e and 7^. 

The agreement between the data and the MC simulation of the leading-proton variables 
is also good for all three samples. As an example, the distributions for the reconstructed 
EPSOFT events as a function of xl, Pti Apipe and A p i ane are compared with those of the 
BPC data in Fig. 3. 
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6 Acceptance 



Figure 4 shows the ranges of px and py accessed by the LPS for six values of xl, using 
the coincidence of any two of the S4, S5 and S6 stations. Here, px and py are the X and 
Y components of the scattered-proton momentum. The region covered is determined by 
the beam-pipe apertures, the shape of the sensitive region of the LPS detectors and the 
magnet strengths; it is limited to xl *S 0.5 and p\ = p\ + p\ ^ 0.7 GeV 2 . Integrated 
over the falling p\ distribution, the LPS geometrical acceptance reaches a maximum for 
xl in the range 0.8 to 0.9. 

The acceptance was computed as the ratio of the number of reconstructed events in a 
bin of a given variable to the number of generated events in that bin. The acceptance 
thus includes the effects of the geometrical acceptance of the apparatus, its efficiency and 
resolution, as well as the event selection and reconstruction efficiencies. Values of the 
acceptance are given in Section 8. 

7 Systematic uncertainties 

The systematic uncertainties were obtained by modifying the requirements and the anal- 
ysis procedures as listed below: 

• the sensitivity to the selection of the proton track was studied by the following proce- 
dure [19]: 

— the track-selection requirements were varied. In particular, the minimum-allowed 
values of A pipe were changed from 0.2 to 0.6 mm and the minimum value of A p i anc 
was varied from 0.1 to 0.3 mm; 

— the positions of some of the elements of the proton beam-line were varied within 
their uncertainties. This is particularly relevant at low xl, where the proton 
momentum is significantly lower than that of the xl ~ 1 protons used in the LPS 
alignment procedure; 

— the LPS detector positions varied slightly from fill to fill in the 94 sample. The 
small deviations of the acceptance implied by these movements were estimated 
by dividing the data into a "low acceptance" and a "high acceptance" sample, 
depending on the positions of the LPS stations; 

• the sensitivity to the remaining selection cuts was also investigated: 

— for the inclusive photoproduction sample, the selection cuts were tightened: the E' e 
range was restricted to 13 < E' e < 16 GeV, corresponding to 195 < W < 215 GeV 
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and the minimum energy deposition in the RCAL was raised to 2 GeV. In addition, 
the correction for the bremsstrahlung background was removed [35]; 

— for the BPC sample, the BPC energy scale was varied by ±1% and the sizes of the 
parameters in the BPC energy calibration were changed within their uncertain- 
ties [24]. The selection limits on the positron-candidate shower width were also 
varied [24]; 

— for the DIS sample, the positron-selection procedure was varied. The cut on the 
scattered-positron energy was changed to 8 GeV and 12 GeV and the size of the 
fiducial region for the impact position of the scattered positron in the SRTD was 
raised by ±0.5 cm in both X and Y; 

— for both the BPC and the DIS samples, the lower limit on 5 was varied between 
32 and 38 GeV and the upper limit between 60 and 68 GeV; the cut on yjs was 
varied by ±0.01; 

— the allowed range of values for the Z coordinate of the vertex was restricted to 
—40 < Z < 40 cm. The effect of removing the vertex requirement was also studied; 

— in addition, for the jet studies, the minimum jet energy was varied between 3.8 
and 4.2 GeV, and the upper limit on rf et was varied between 2 and 2.4. 

The total systematic uncertainty on the cross sections, obtained by summing all the above 
contributions in quadrature, totalled about ±20% at xl ~ 0.65, decreasing to ±(10- 
15)% for xl ^ 0.75. The dominant contributions are those related to the track-selection 
requirements in the LPS. 



8 The ratio method 

In the following, several results are presented in terms of the ratios, r LP ^ and r LP ^ 3 \ of 
the cross section for production of leading protons to the cross section for inclusive e + p 
scattering; these ratios are evaluated in bins of x and Q 2 (r LP ( 2 )), or in bins of x, Q 2 and Xl 
(r Lp (3)). They are obtained from the measured fraction of the events, in a given bin, that 
have a leading proton, N hP /N. In this fraction, the acceptance corrections related to the 
positron selection procedure cancel, and so do the corresponding systematic uncertainties. 
The only remaining correction to apply is that for the LPS acceptance, eLPs- 

The ratio r LP ^ was thus obtained as 

r LP W(x,Q 2 ) ArLP ^ C ? 2 ) 1 



N(x,Q 2 ) e LPS 



Averaged over the region 0.6 < xl < 0.97 and y\ < 0.5 GeV 2 , e^ps is approximately 23%; 
over the region 0.6 < xl < 0.97 and p T < 0.04 GeV 2 , eip$ ~ 51%. 
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From Eqs. (1) and (2), it is apparent that the cross section ratio r LP ( 2 ) is also equal to 
the ratio of the proton tagged and inclusive structure functions: 

F hp W(r O 2 ) 

^,0) = *L-gQ), (3 ) 
where F 2 LP( - 2 ^ is obtained from by integration over the measured xl and p\ ranges: 

PTmax 0.97 

Lp ( 2 )^ rfi\ _ f jji f p lp ( 4 )^ n2 ^2 



F™(*,Q 2 ) = J dp 2 T J dx L ( >(x,Q 2 ,x l ,Pt)- 

0.6 

The radiative corrections and the contributions from F L are assumed to be the same for 
the inclusive and the proton-tagged reactions. 

The ratio r LP ( 3 \x, Q 2 , x L ) is defined in analogy to r LP ( 2 \x, Q 2 ): 



^\x,Q\x L ) = Nh ^ Q ^ L) 



N(x,Q 2 ) e LP5 (x i )Ax i ' 
where Axl indicates the size of the xl bins. In analogy with Eq. (3), 



(4) 



r ^,Q^ l) = ggg|p) , (5 ) 

where F^ P( - 3 ' ) (a;, Q 2 , Xl) differs from F^ (x, Q 2 ) in that no integration over xl is per- 
formed. 

The ratios r LF ^ and r LF ^ can also be interpreted in terms of the virtual photon-proton 
cross section for the process 7*p — > Xp and the total virtual photon-proton cross section, 
cr t ot- For example, the ratio r LF ^ can be written as 

Pt max ^ 

LP(3)/ n 2 \ _ _ f i 2 ^ C-y*p^Xp _ _ _dfj*p-^Xp_ 



r^(x,Q 2 ,x L ) = — / dp T ——- 2 

0"tot J dxidptp cr to t dx L 

o 

where the virtual photon-proton cross section, d 2 a 1 * p ^x P /dxLdp 2 ^, is related to the positron- 
proton cross section, d 4 a e + p ^ e +x p /dQ 2 dxdxLdp 2 ^, by 

d 4 a e + p ^ e +xp - r d 2 a^* p ^xp 



dQ 2 dxdx^dp\ dx^dp\ ' 
where T = {a/xQ 2 rc)[l + (1 — y 2 )], is the photon flux factor and a is the fine-structure 



constant. 
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9 Models 



The data were tested against the hypothesis of vertex factorisation, a very general feature 
of hadron-hadron interactions [1] which relates reactions with different beam particles to 
their respective total cross sections. In particular, in the reaction 7*p — > Xp, the 7*-X 
and p-p vertices factorise if the amplitude for the reaction can be written as the product 
of two vertex functions, G^xix, Q 2 ) and G pp (xl,Pt). In this case, the cross section as a 
function of the lepton variables x and Q 2 should be independent of the baryon variables 
Xl and p\, and vice versa. 

The data were also compared to the following specific models: 

• the LUND string-fragmentation model as implemented in JETSET [47] and used in 
DJANGO [48], in which leading baryons originate from the hadronisation of the target; 

• the soft-colour-interaction model (SCI) [17], as implemented in LEPTO 6.5 [49], where 
leading protons are obtained from standard DIS events by means of a non-perturbative 
redistribution of colour among the fragmenting partons; 

• the Regge model of Szczurek et al. [16], which assumes a superposition of Pomeron, 
Reggeon and pion exchanges; 

• the QCD-inspired model of Duraes et al. [50], developed in the framework of the 
interacting-gluon model [51], which assumes that high-energy hadron-hadron collisions 
are dominated by multiple incoherent gluon-gluon interactions. The valence quarks 
that do not take part in the interaction give rise to the leading baryons. The extension 
to ep collisions is made in the VDM framework. 

10 Results 

10.1 Leading proton events with a forward large rapidity gap 

Some indication of the production mechanism of leading protons can be obtained from 
the rapidity distribution of the hadronic final-state particles. In particular, events of 
diffractive origin, i.e. due to Pomeron exchange, are characterised by a gap in the rapidity 
distribution in the forward direction. 

Figure 5a) shows the distribution of the DIS events in the (^ max , xl) plane, where ?7 max 
is the pseudorapidity of the most-forward energy deposit of at least 400 MeV in the 
CAL. The accumulation of events at xl ~ 1, which mostly have r] max < 2.5, is due 
to diffractive events [52], e + p —> e + Xp, in which the final-state proton remains intact 
and carries approximately the same momentum as the incoming proton. Events with 
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?7max < 2.5 and Xl ^ 0.97 are ascribed to double diffractive dissociation, e + p — > e + XN, 
where the proton dissociates into the state N, with mass Mjv. Although N is produced 
at xl ~ 1, the proton from the decay of N has a lower value of xl- When both M N and 
Mx are small, the systems X and N are separated by a large gap in pseudorapidity. If 
M N < 4-5 GeV, only the proton from the system N is measured, while the other particles 
escape undetected down the beam pipe. The topology of a doubly dissociative event is 
thus characterised by a rapidity gap in the forward region in conjunction with a low-x^ 
proton. 

To select a sample of diffractive events, the requirement ?7 max < 2.5 was imposed. Fig- 
ure 5b) shows the fraction of BPC and DIS events with a leading proton that also have 
Vm&x < 2.5 (see also Table 1). Events with a large rapidity gap dominate for xl ~ 1. 
For 0.6 < xl < 0.97, the fraction of leading proton events with a large rapidity gap is 
less than 10% in any given bin, and is only weakly dependent on Q 2 and xl- This result, 
which indicates that diffraction is not the main mechanism responsible for production 
of leading protons in this region, is consistent with the predictions of the Regge-based 
models [16, 53]. 

10.2 Momentum spectra of leading protons 
10.2.1 Longitudinal-momentum spectra 

The normalised cross-section r LP ^ = (l/a tot )-da^* p ^xp/dxL for the reaction e + p — > e + Xp 
with a leading proton having xl > 0.6 and p\ < 0.5 GeV 2 is shown in Fig. 6 and 
given in Table 2 for the BPC sample, integrated over the range 0.1 < Q 2 < 0.74 GeV 2 , 
85 < W < 258 GeV, 1.5 x 10~ 6 < x < 1.0 x 10~ 4 , and for the DIS sample, integrated 
over the region 3 < Q 2 < 254 GeV 2 , 45 < W < 225 GeV, 1.2 x 10~ 4 < x < 4 x 10~ 2 . 
These results are compared with those from the reaction pp — > pX at y^s = 19.6 GeV [54], 
integrated over the same p\ region and normalised to the corresponding inelastic cross 
section. For xl ^ 0.9, the fraction of events with a leading proton is consistent for the 
pp and 7*p data sets, in accord with vertex factorisation. A dependence on the centre- 
of-mass energy is apparent for xl > 0.9, as expected from Regge parametrisations of the 
pp — > pX data [55, 56]. 

Figure 7 and Table 3 present the photoproduction, BPC and DIS data for the lower p\ 
region, p\ < 0.04 GeV 2 and 0.6 < xl < 0.95, where the upper cut on xl, which removes 
the diffractive events, is set by the LPS acceptance for this p\ range. The fraction of 
events with a leading proton is approximately the same in all three regimes. The pp 
data [54] for p\ < 0.05 GeV 2 again agree with the ep data for xl ^ 0.9. The present 
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photoproduction results, however, are significantly higher than those found by HI [7] in 
similar ranges of Q 2 , W and p\. 

Figure 8 compares the DIS data to the specific models of Section 9. The standard DIS 
Monte Carlo generator DJANGO [48] predicts a stronger decrease of the cross section 
with xl than that observed and, in addition, has no diffractive peak. It also substantially 
underestimates the rate of leading proton production in the measured xl range. The SCI 
model [17], as implemented in LEPT06.5 [49], also falls below the data, even though it 
generates a larger number of leading protons than DJANGO and has a peak at xl = 1. 
These two models are ruled out by the data. 

The QCD-inspired model of Duraes et al. [50, 57] is in better agreement with the data, 
but is too low in the diffractive peak region. Nevertheless, the similarity to the data is 
remarkable, given the small number of free parameters in the model. 

The Regge-based calculation of Szczurek et al. [16] agrees with the data, although it is 
somewhat too low at small values of xl- In this approach, leading proton production 
for 0.6 < xl < 0.9 is dominated by isoscalar Reggeon exchange; diffractive processes, 
due to Pomeron exchange, become increasingly important as II approaches unity. The 
contribution of pion exchange, including the production of A baryons, is also evaluated. 
The normalisation of the Reggeon contribution has a large theoretical uncertainty [16,58] 
and the present data suggest that this contribution should be increased with respect to 
the assumption made by Szczurek et al. [16]. 

10.2.2 Transverse-momentum spectra 

The cross-sections da^p^xp/ dp\ are shown in Fig. 9 for the BPC sample integrated over 
the range 0.1 < Q 2 < 0.74 GeV 2 , 85 < W < 258 GeV, 1.5 x 10~ 6 < x < 1.0 X 10~ 4 
for different x L selections. Similar distributions for the DIS sample are shown in Fig. 10, 
integrated over the region 3 < Q 2 < 254 GeV 2 , 45 < W < 225 GeV, 1.2 x 10~ 4 < x < 
4 x 10~ 2 . In all regions, the form da^* v ^xpl dp\ = Ae~ bp T fits satisfactorily the data. 
The values of the slope-parameters b obtained for each xl bin are, within uncertainties, 
independent of shown in Fig. 11 (see also Table 4). 

The BPC and DIS data together indicate that b is independent of Q 2 and xl- The 
mean value of b for the BPC data is (b) = 6.6 ± 0.6 (stat.) ± 0.8 (syst.) GeV -2 and 
(b) = 6.9 ± 0.2 (stat.) ± 0.8 (syst.) GeV~ 2 for the DIS data for 0.6 < x L < 0.97. Also 
plotted in Fig. 11a) is the result obtained for diffractive photoproduction [59], which 
is consistent with the values found at higher Q 2 . In addition, the present results are 
compatible with the pp — > pX data [54], also shown in Fig. 11a). This, together with 
the fact that b is approximately Q 2 -independent, provides additional support for vertex 
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factorisation. 

The predictions of Szczurek et al. [16] are in accord with the transverse-momentum data, 
as shown in Fig. lib). In this model, no Q 2 dependence is expected for b [60]. The 
DJANGO program (not shown) has an effective slope b pa 4 GeV" 2 . LEPT06.5 (also 
not shown) has b ~ 3.5 GeV -2 . Both are independent of xl and significantly below the 
values measured. The model of Duraes et al. [50] does not make explicit predictions for 
the transverse-momentum distribution. 

10.3 The proton-tagged structure-function 

The acceptance-corrected fraction of events with a leading proton is used to measure the 
proton-tagged structure- function -F 2 LP , as discussed in Section 8. To select a predominantly 
non-diffractive sample, the cut xl < 0.97 was imposed. 

Figures 12 and 13 show r LP ® (x, Q 2 , x^), determined using Eq. (4), for the BPC and DIS 
samples, respectively, in several x and Q 2 bins for p\ < 0.5 GeV 2 . The data are also given 
in Tables 5-7. Only a weak dependence on x and Q 2 is apparent, indicating that F 2 LP ® 
has approximately the same x and Q 2 dependence as F 2 (x, Q 2 ). The data exhibit a weak 
xl dependence, as already seen in Fig. 6. Figure 14 shows r LP ^ for fixed Q 2 values as a 
function of x; the results again have little x dependence (see also Table 8). 

Figure 15a) and Table 9 present the BPC and DIS data for 0.6 < xl < 0.97 and p\ < 
0.5 GeV 2 , averaged over x for different Q 2 ranges, (r LP( - 2 \Q 2 )) . The leading proton 
yield increases by approximately 20%, from (r LP ( 2 ^) pa 0.12 to (r LP ( 2 )) pa 0.15, when Q 2 
varies from ~ 0.25 GeV 2 (the average value of Q 2 for the BPC sample) to 100 GeV 2 , 
indicating a modest but definite breakdown of vertex factorisation. Figure 15b) and 
Table 10 present the points of Fig.l5a) normalised to (r hp ( 2 \Q 2 = 0.25 GeV 2 )). The 
results for the restricted p\ range, p\ < 0.04 GeV 2 , are also shown. The breaking of vertex 
factorisation is approximately the same for p\ < 0.5 GeV 2 and for p\ < 0.04 GeV 2 . An 
effect of similar size was measured for leading neutron production [8]; the corresponding 
data, normalised to the value at Q 2 = 0.25 GeV 2 , are also shown in Fig. 15b). The 
neutron data are measured for scattering angles less than 0.8 mrad, corresponding to 
p% < 0.43 • x\ GeV 2 . 

This Q 2 dependence of the proton yield can be qualitatively ascribed to absorptive effects 
in the 7*p system [11]. The transverse size of the virtual photon decreases with increasing 
Q 2 , reducing the likelihood that the produced baryon rescatters on the hadronic compo- 
nent of the virtual photon. 

The data of Fig. 14 are presented in Fig. 16 in terms of F 2 , obtained by multiplying 
r Lp (2) j-^ j? 2 _ p or the BPC region, a parameterisation of the ZEUS F 2 results [25] was 
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used. For the DIS region, the parameterisation of Botje [61] was used. Since r LP ( 2 ) is 
approximately independent of Q 2 and x, is approximately proportional to F%. As 

indicated in the figure, the proportionality constant between and F 2 is (r LP ( 2 )) ps 

0.13. 

The HI collaboration has published [5] a study of leading proton production for p\ < 
0.04 GeV 2 , 0.7 < x L < 0.9, 2 < Q 2 < 50 GeV 2 and 6 x 10" 5 < x < 6 x 10~ 3 . The present 
analysis was repeated in the region of overlap with the HI data set. A comparison of the 
HI and ZEUS results is presented in Fig. 17. The agreement is good. 

10.4 Leading protons with associated dijet production 

The ratio r^p of the yield of leading proton DIS events with associated dijet production 
to the inclusive yield of leading proton events is presented as a function of xl and p\ 
in Fig. 18 and Tables 13-14. The LPS acceptance, as well as the scattered positron 
acceptance, cancels in this ratio. The data are shown for the range 0.6 < xl < 0.97 and 
p\ < 0.5 GeV 2 . No significant deviation from a flat behaviour is seen as a function of xl, 
although there is some p\ dependence. The results of this exploratory study thus suggest 
that the longitudinal- and transverse-momentum distributions of the leading proton are 
largely insensitive to the presence of a second hard scale, given by the transverse energy 
of the jets. 

Figure 19 and Tables 15-17 present the fraction of the dijet events with a leading proton, 
r^ P , plotted as a function of Et, x and Q 2 . In this case, all corrections cancel with the 
exception of that due to the LPS acceptance, which is, however, independent of E T , Q 2 
and x. The ratio is approximately independent of these variables and its value is consistent 
with that of r LP ( 2 \ This suggests that the Et, Q 2 and x dependences of the dijet cross 
section are unaffected by the requirement of a leading proton, and that the fraction of 
dijet events with a leading proton is the same as the fraction of inclusive events with a 
leading proton. 

10.5 Summary 

Events of the type e + p — > e + Xp with a final-state proton with xl > 0.6 have been studied 
in e + p collisions at HERA. The analyses used a photoproduction sample (Q 2 < 0.02 GeV 2 ), 
a low-Q 2 sample (0.1 < Q 2 < 0.74 GeV 2 ) and a DIS sample (3 < Q 2 < 254 GeV 2 ). 

For events with a leading proton in the range 0.6 < xl < 0.97 and p\ < 0.5 GeV 2 , the 
main features of the data can be summarised as follows: 
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• less than 10% of the leading proton events in any given Xl bin exhibit a large rapidity 
gap (r/ max < 2.5), indicating that diffraction is not the main mechanism responsible 
for the production of leading protons in this region; 

• the proton xl spectrum is only a weak function of xl for £l^0.97; 

• the p\ dependence of the cross section is well described by an exponential function, 
with a slope approximately independent of xl and Q 2 , b ~ 7 GeV~ 2 . The slope is also 
consistent with the value measured for pp collisions; 

• the x and Q 2 dependence of the semi-inclusive structure function, F 2 LP , is similar to 
that of F 2 , independently of Xl- However, F 2 LP grows with Q 2 slightly faster than F 2 , 
resulting in a yield of leading protons about 20% larger at Q 2 = 100 GeV 2 than at 
Q 2 ~ 0. A similar effect was observed for leading neutron production; 

• the shapes of the Xl and p\. spectra are largely unaffected by requiring two jets within 
the hadronic final state X. The Q 2 , x and Et dependences of the dijet cross section 
are also broadly consistent for leading proton events and inclusive events. 

The main features of the data are reproduced by a Regge model assuming a superposition 
of Pomeron, Reggeon and pion trajectories. The D JANGO and SCI models are ruled out 
by the data. 

For 0.6 < Xl < 0.9, the proton spectrum for (virtual-) photon-proton collisions is con- 
sistent with the results found in proton-proton reactions at lower centre-of-mass energy. 
The fraction of the events with a leading proton is approximately the same for the 7*p 
and pp data, in agreement with vertex factorisation. 

In the xl region explored, a modest violation of vertex factorisation is observed. Never- 
theless, the results of this paper indicate that the properties of the final-state proton are 
largely independent of those of the virtual photon. 
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Table 1: Fraction of events with rj max < 2.5. Statistical uncertainties are given; 
systematic uncertainties mostly cancel in the ratio. 



XL 


^LP^max 


< 2.5)/iY LP 


BPC 


0.63 


0.0o4 


± 





030 


0.66 


0.116 


± 





033 


0.69 


0.0 10 


± 


() 


mo 
018 


0.72 


0.070 


± 





020 


0.75 


0.083 


± 





021 


0.78 


0.099 


± 





020 


0.81 


0.077 


± 





018 


0.84 


0.044 


± 





014 


0.87 


0.066 


± 





017 


0.90 


0.042 


± 





014 


0.93 


0.011 


± 





008 


0.96 


0.040 


± 





016 


0.99 


0.355 


± 





064 


1.00 


0.812 


± 





036 


DIS 


0.63 


0.0674 


± 





0093 


0.66 


0.0567 


± 





0081 


0.69 


0.0728 


± 





0072 


0.72 


0.0603 


± 





0063 


0.75 


0.0543 


± 





0057 


0.78 


0.0507 


± 





0050 


0.81 


0.0597 


± 





0047 


0.84 


0.0495 


± 





0049 


0.87 


0.0364 


± 





0040 


0.90 


0.0421 


± 





0047 


0.93 


0.0304 


± 





0053 


0.96 


0.0805 


± 





0132 


0.99 


0.3445 


± 





0231 


1.00 


0.7683 


± 





0164 
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Table 2: The normalised cross-section (l/cr tot ) ■ da y * p ^,xp/dxL for the BPC and 
DIS data in the region p\ < 0.5 GeV 2 . The two rightmost values indicate the 
statistical and systematic uncertainties, respectively. 



XL 



(l/otot) " da^p^xp/dxi 



BPC 



0.65 
0.75 
0.83 
0.93 
0.99 



0.339 ± 0.022 ± 0.050 
0.399 ± 0.019 ± 0.048 
0.330 ± 0.014 ± 0.049 
0.331 ± 0.021 ± 0.040 
3.60 ± 0.37 ± 0.54 



DIS 


0.62 


0.413 


± 





019 


± 0.091 


0.65 


0.431 


± 





017 


± 0.065 


0.69 


0.462 


± 





016 


± 0.056 


0.71 


0.445 


± 





015 


± 0.066 


0.75 


0.421 


± 





012 


± 0.088 


0.77 


0.433 


± 





012 


± 0.048 


0.81 


0.379 


± 





010 


± 0.046 


0.83 


0.359 


± 





009 


± 0.029 


0.87 


0.368 


± 





010 


± 0.037 


0.89 


0.333 


± 





010 


± 0.030 


0.93 


0.289 


± 





012 


± 0.026 


0.95 


0.46 ± 0.03 ± 0.11 


0.99 


2.48 


± 0.12 ± 0.37 
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Table 3: The normalised cross-section (l/cr tot ) ■ da^p^xp/dxi for the photopro- 
duction, BPC and DIS data for p\ < 0.04 GeV 2 . The two rightmost values indicate 
the statistical and systematic uncertainties, respectively. 



XL 


(l/o"tot) • da r * p ^x P /dx L 


Photoproduction 


0.64 


0.110 ± 0.017 ± 0.022 


0.70 


0.081 ± 0.008 ± 0.016 


0.76 


0.079 ± 0.006 ± 0.012 


0.82 


0.090 ± 0.006 ± 0.012 


0.88 


0.080 ± 0.006 ± 0.012 


BPC 


0.65 


0.084 ± 0.007 ± 0.013 


0.75 


0.099 ± 0.005 ± 0.012 


0.83 


0.081 ± 0.005 ± 0.012 


0.93 


0.0697 ± 0.0100 ± 0.0084 


DIS 


0.62 


0.111 ± 0.008 ± 0.024 


0.65 


0.105 ± 0.006 ± 0.016 


0.69 


0.110 ± 0.005 ± 0.013 


0.71 


0.103 ± 0.004 ± 0.015 


0.75 


0.095 ± 0.003 ± 0.020 


0.77 


0.102 ± 0.003 ± 0.011 


0.81 


0.092 ± 0.003 ± 0.011 


0.83 


0.0910 ± 0.0030 ± 0.0073 


0.87 


0.0910 ± 0.0040 ± 0.0091 


0.89 


0.0820 ± 0.0050 ± 0.0075 
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Table 4: The slopes, b, from fits of the functional form e~ bPr to da^p^xp/dPr f or 
leading protons as a function of xl for the BPC and DIS data samples. The two 
rightmost values indicate the statistical and systematic uncertainties, respectively. 



XL 


b (GeV" 


- 2 ) 


BPC 


0.65 


7.5 


± 1.2 


+0.1 
-1.9 


0.75 


7.6 ± 1.0 


+1.2 
-1.1 


0.83 


6.2 ± 1.2 


+1.2 
-0.7 


0.93 


4.5 ± 1.2 


+1.5 
-0.8 


DIS 


0.62 


7.4 


± 


1.2 


+1.0 
-1.3 


0.65 


7.7 


± 


0.8 


+2.0 
-1.2 


0.69 


7.3 


± 


0.6 


+0.8 
-0.8 


0.71 


6.8 


± 


0.6 


+0.3 
-0.6 


0.75 


6.1 


± 


0.6 


+1.2 
-1.2 


0.77 


6.8 


± 


0.7 


+0.4 
-0.4 


0.81 


7.9 


± 


0.9 


+0.8 
-0.9 


0.83 


7.3 


± 


0.8 


+0.6 
-0.9 


0.87 


7.2 


± 


0.7 


+ 1.1 
-1.8 


0.89 


6.6 


± 


0.7 


+0.6 
-0.5 


0.93 


5.9 


± 


0.9 


+1.2 
-1.0 


0.95 


4.4 


± 


1.7 


+1.5 
-0.4 


0.99 


7.0 


± 


0.9 


+1.9 
-1.1 
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Table 5: The ratio r hP ^ = F 2 LP(3) /-P 2 as a function of x L , x and Q 2 (BPC 
sample), for protons with p\ < 0.5 GeV 2 . The two rightmost values indicate the 
statistical and systematic uncertainties, respectively. 



Q 2 (GeV 2 ) 


X 


XL 


r LP(3) 


0.2 


7.4E-06 


o 


67 


0.310 


± 


o 


050 


± 


0.056 


0.2 


7.4E-06 


o 


79 


0.408 


± 


o 


042 


± 


0.049 


0.2 


7.4E-06 


o 


91 


0.367 


± 


o 


049 


± 


0.037 


0.2 


4.9E-06 





67 


0.308 


± 


o 


052 


± 


0.055 


0.2 


4.9E-06 


o 


79 


0.331 


± 


o 


039 


± 


0.040 


0.2 


4.9E-06 





91 


0.361 


± 


o 


050 


± 


0.036 


0.2 


3.5E-06 


o 


67 


0.281 


± 


o 


063 


± 


0.051 


0.2 


3.5E-06 


o 


79 


0.333 


± 


o 


049 


± 


0.040 


0.2 


3.5E-06 


o 


91 


0.377 


± 


o 


065 


± 


0.038 


0.4 


2.6E-05 





67 


0.357 


± 





031 


± 


0.064 


0.4 


2.6E-05 





79 


0.356 


± 





022 


± 


0.043 


0.4 


2.6E-05 





91 


0.317 


± 





026 


± 


0.032 


0.4 


1.3E-05 





67 


0.356 


± 





061 


± 


0.064 


0.4 


1.3E-05 





79 


0.405 


± 





047 


± 


0.049 


0.4 


1.3E-05 





91 


0.265 


± 





047 


± 


0.027 


0.4 


8.8E-06 





67 


0.357 


± 





084 


± 


0.064 


0.4 


8.8E-06 





79 


0.317 


± 





057 


± 


0.038 


0.4 


8.8E-06 





91 


0.298 


± 





068 


± 


0.030 


0.6 


4.3E-05 





67 


0.438 


± 





063 


± 


0.079 


0.6 


4.3E-05 





79 


0.386 


± 





043 


± 


0.046 


0.6 


4.3E-05 





91 


0.354 


± 





051 


± 


0.035 
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Table 6: The ratio r LP ^ = F 2 LP< - 3 ^/F 2 as a function ofxz, x andQ 2 (DIS sample, 
up to Q 2 = 12 GeV 2 ), for protons with p\ < 0.5 GeV 2 . The two rightmost values 
indicate the statistical and systematic uncertainties, respectively. 



Q 2 (GeV 2 ) 


X 


%L 


r LP(3) 


4.0 


1 


5E 


04 





67 


0.344 


± 0.033 


± 0.062 


4.0 


1 


5E 


04 





79 


0.345 


± 0.023 


± 0.041 


4.0 


1 


5E 


04 





91 


0.326 


± 0.029 


± 0.033 


4.0 


2 


5E 


04 





67 


0.327 


± 0.0.35 


± 0.059 


4.0 


2 


5E 


04 





79 


0.353 


± 0.026 


± 0.042 


4.0 


2 


5E 


04 





91 


0.313 


± 0.031 


±0.031 


4.0 


4 


-IE 


04 





67 


0.367 


± 0.035 


± 0.066 


4.0 


4 


4E 


04 





79 


0.335 


± 0.024 


± 0.040 


4.0 


4 


4E 


04 





91 


0.305 


± 0.029 


±0.030 


4.0 


9 


8E 


04 





67 


0.443 


± 


0.040 


± 0.080 


4.0 


9 


8E 


04 





79 


0.354 


± 


0.026 


± 0.042 


4.0 


9 


8E 


04 





91 


0.303 


± 


0.030 


± 0.030 


8.0 


2 


5E 


04 





67 


0.398 


± 


0.031 


± 0.072 


8.0 


2 


5E 


04 





79 


0.359 


± 


0.021 


± 0.043 


8.0 


2 


5E 


04 





91 


0.382 


± 


0.028 


± 0.038 


8.0 


4 


4E 


04 





67 


0.387 


± 


0.028 


± 0.070 


8.0 


4 


4E 


04 





79 


0.409 


± 


0.020 


± 0.049 


8.0 


4 


4E 


04 





91 


0.335 


± 


0.023 


± 0.033 


8.0 


9 


8E 


04 





67 


0.451 


± 


0.031 


± 0.081 


8.0 


9 


8E 


04 





79 


0.420 


± 


0.021 


± 0.050 


8.0 


9 


8E 


04 





91 


0.375 


± 


0.026 


± 0.038 


8.0 


2 


2E 


03 





67 


0.515 


± 


0.050 


± 0.093 


8.0 


2 


2E 


03 





79 


0.380 


± 


0.031 


± 0.046 


8.0 


2 


2E 


03 





91 


0.421 


± 


0.041 


± 0.042 


12.0 


4 


4E 


04 





67 


0.355± 


0.029 ± 0.064 


12.0 


4 


IE 


04 





79 


0.379 


± 


0.022 ± 0.046 


12.0 


4 


IE 


04 





91 


0.396 


± 


0.028 ± 0.040 


12.0 


9 


8E 


04 





67 


0.416 


± 


0.029 ± 0.075 


12.0 


9 


8E 


04 





79 


0.387 


± 


0.020 ± 0.046 


12.0 


9 


8E 


04 





91 


0.369 


± 


0.025 ± 0.037 


12.0 


2 


2E 


03 





67 


0.437 


± 


0.039 ± 0.079 


12.0 


2 


2E 


03 





79 


0.407 


± 


0.027 ± 0.049 


12.0 


2 


2E 


03 





91 


0.445 


± 


0.0.36 ± 0.044 


12.0 


4 


0E 


03 





67 


0.462 


± 


0.067 ± 0.083 


12.0 


4 


OE 


03 





79 


0.427 


± 


0.046 ± 0.051 


12.0 


4 


OE 


03 





91 


0.349 


± 


0.053 ± 0.035 
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Table 7: The ratio r LP ^ = F^ F /F 2 as a function ofxi, x and Q 2 (DIS sample, 
for Q 2 > 12 GeV 2 ;, for prot ons with p\ < 0.5 GeV 2 . The two rightmost values 
indicate the statistical and systematic uncertainties, respectively. 



Q 2 (GeV ) 


X 


XL 


r LP(3) 


21.0 


9 


8E- 


04 





67 





423 


± 





027 


± 


0.076 


21.0 


9 


8E- 


04 





79 





387 


± 





019 


± 


0.046 


21.0 


9 


8E- 


04 





91 





360 


± 





023 


± 


0.036 


21.0 


2 


2E- 


03 





67 





437 


± 





034 


± 


0.079 


21.0 


2 


2E- 


03 





79 





427 


± 





024 


± 


0.051 


21.0 


2 


2E- 


03 





91 





379 


± 





029 


± 


0.038 


21.0 


4 


0E- 


03 





67 





486 


± 





049 


± 


0.087 


21.0 


4 


OE- 


03 





79 





390 


± 





031 


± 


0.047 


21.0 


4 


OE- 


03 





91 





361 


± 





038 


± 


0.036 


21.0 


5 


9E- 


03 





67 





474 


± 





066 


± 


0.085 


21.0 


5 


9E 


03 





79 





372 


± 





042 


± 


0.045 


21.0 


5 


9E- 


03 





91 





466 


± 





059 


± 


0.047 


46.0 


2 


2E- 


03 





67 





452 


± 





043 


± 


0.081 


46.0 


2 


2E- 


03 





79 





372 


± 





028 


± 


0.045 


46.0 


2 


2E- 


03 





91 





490 


± 





040 


± 


0.049 


46.0 


4 


OE- 


03 





67 





444 


± 





055 


± 


0.080 


46.0 


4 


OE- 


03 





79 





433 


± 





039 


± 


0.052 


46.0 


4 


OE- 


03 





91 





360 


± 





045 


± 


0.036 


46.0 


5 


9E 


03 





67 





435 


± 





064 


± 


0.078 


46.0 


5 


9E- 


03 





79 





348 


± 





041 


± 


0.042 


46.0 


5 


9E- 


03 





91 


0.436 


± 


0.059 


±0.044 


46.0 


1 


OE- 


02 





67 





453 


± 





050 


± 


0.081 


46.0 


1 


OE- 


02 





79 





412 


± 





034 


± 


0.049 


46.0 


1 


OE- 


02 





91 





409 


± 





043 


± 


0.041 


130.0 


5 


9E 


03 





67 





441 


± 





094 


± 


0.079 


130.0 


5 


9E- 


03 





79 





327 


± 





058 


± 


0.039 


130.0 


5 


9E 


03 





91 





528 


± 





094 


± 


0.053 


130.0 


1 


OE- 


02 





67 





408 


± 





060 


± 


0.073 


130.0 


1 


OE- 


02 





79 





505 


± 





048 


± 


0.061 


130.0 


1 


OE- 


02 





91 





358 


± 





051 


± 


0.036 


130.0 


2 


5E- 


02 





67 





421 


± 





069 


± 


0.076 


130.0 


2 


5E- 


02 





79 





430 


± 





050 


± 


0.052 


130.0 


2 


5E- 


02 





91 





403 


± 





062 


± 


0.040 
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Table 8: The ratio r LP ^ = F 2 LP(2) /F 2 as a function of x for fixed Q 2 values, for 
protons with 0.6 < xl < 0.97 and p\ < 0.5 GeV 2 . The statistical uncertainty is 
given. A fully correlated systematic uncertainty of ±13% is not included. 



Q 2 (GeV ) 


X 


r LP(2) 


0.20 


3 


5E- 


06 


0.118 


± 





012 


0.20 


4 


9E- 


06 


0.1177 


± 





0093 


0.20 


7 


4E- 


06 


0.1309 


± 





0095 


0.36 


8 


8E- 


06 


0.113 


± 





014 


0.36 


1 


3E 


05 


0.124 


± 





010 


0.36 


2 


6E- 


05 


0.1211 


± 





0092 


0.60 


4 


3E- 


05 


0.137 


± 





010 


4.0 


9 


8E- 


04 


0.1281 


± 





0063 


4.0 


4 


4E- 


04 


0.1186 


± 





0058 


4.0 


2 


5E- 


04 


0.1190 


± 





0062 


4.0 


1 


5E- 


04 


0.1206 


± 





0056 


8.0 


2 


2E- 


03 


0.1507 


± 





0079 


8.0 


9 


8E- 


04 


0.1473 


± 





0052 


8.0 


4 


4E- 


04 


0.1358 


± 





0048 


8.0 


2 


5E- 


04 


0.1332 


± 





0053 


12.0 


4 


0E- 


03 


0.147 


± 





Oil 


12.0 


2 


2E- 


03 


0.1511 


± 





0067 


12.0 


9 


8E- 


04 


0.1382 


± 





0049 


12.0 


4 


4E- 


04 


0.1345 


± 





0053 


21.0 


5 


9E- 


03 


0.151 


± 





Oil 


21.0 


4 


OE- 


03 


0.1438 


± 





0077 


21.0 


2 


2E- 


03 


0.1477 


± 





0058 


21.0 


9 


8E- 


04 


0.1377 


± 





0046 


46.0 


1 


OE- 


02 


0.1496 


± 





0084 


46.0 


5 


9E 


03 


0.140 


± 





Oil 


46.0 


4 


OE- 


03 


0.1472 


± 





0092 


46.0 


2 


2E- 


03 


0.1512 


± 





0072 


130.0 


2 


5E- 


02 


0.149 


± 





012 


130.0 


1 


OE- 


02 


0.156 


± 





Oil 


130.0 


5 


9E- 


03 


0.147 


± 





016 
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Table 9: The average ratio (r LP ( 2 )) = F 2 LP(2) /F 2 as a function ofQ 2 for 0.6 < Xl < 
0.97 and Pt < 0.5 GeV 2 . The statistical uncertainty is given. A fully correlated 
systematic uncertainty of 13% is not included. 



Q 2 (GeV 2 ) 


( r LP(2)} 


0.29 


0.1230 ± 0.0033 


5.21 


0.1312 ± 0.0020 


16.60 


0.1397 ± 0.0021 


69.00 


0.1471 ± 0.0036 



Table 10: The average ratio (r LP ( 2 )) as a function of Q 2 for two different p\ 
ranges normalised to the value at Q 2 = 0.25 GeV 2 . The statistical uncertainty is 
given; systematic errors mostly cancel in the ratio. 



Q 2 (GeV 2 ) 


(r LP M(Q 2 ))/{r LP W(Q 2 = 0.25 GeV 2 )) 




p 2 T < 0.04 GeV 2 


p 2 T < 0.5 GeV 2 


0.002 


0.941 ± 0.033 




0.29 


1.000 ± 0.000 


1.000 ± 0.000 


5.21 


1.062 ± 0.022 


1.067 ± 0.031 


16.60 


1.115 ± 0.023 


1.136 ± 0.031 


69.00 


1.152 ± 0.039 


1.196 ± 0.037 
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Table 11: The structure-function F 2 as a function of x for 0.6 < Xl < 
0.97 and p\ < 0.5 GeV 2 . The statistical uncertainty is given. A fully correlated 
systematic uncertainty of ±13% is not included, nor is the uncertainty of the F 2 
parametrisations used. 



X 


Q (GeV ) 


Fr 


LP(2) 




3 


5E- 


06 


0.2 


0.0286 


± 





0029 


4 


9E- 


06 


0.2 


0.0277 


± 





0022 


7 


4E- 


06 


0.2 


0.0296 


± 





0021 


8 


8E- 


06 


0.4 


0.0391 


± 





0048 


1 


3E- 


05 


0.4 


0.0412 


± 





0035 


2 


6E 


05 


0.4 


0.0378 


± 





0029 


4 


3E 


05 


0.6 


0.0559 


± 





0042 


1 


5E- 


04 


4.0 


0.1126 


± 





0052 


2 


5E- 


04 


4.0 


0.0991 


± 





0052 


4 


4E- 


04 


4.0 


0.0872 


± 





0043 


9 


8E- 


04 


4.0 


0.0793 


± 





0039 


2 


5E- 


04 


8.0 


0.1518 


± 





0060 


4 


4E- 


04 


8.0 


0.1347 


± 





0048 


9 


8E- 


04 


8.0 


0.1201 


± 





0042 


2 


2E- 


03 


8.0 


0.1014 


± 





0053 


4 


4E- 


04 


12.0 


0.1563 


± 





0062 


9 


8E- 


04 


12.0 


0.1303 


± 





0046 


2 


2E- 


03 


12.0 


0.1157 


± 





0051 


4 


0E- 


03 


12.0 


0.0969 


± 





0073 


9 


8E- 


04 


21.0 


0.1552 


± 





0052 


2 


2E- 


03 


21.0 


0.1325 


± 





0052 


4 


OE- 


03 


21.0 


0.1093 


± 





0059 


5 


9E- 


03 


21.0 


0.1029 


± 





0074 


2 


2E- 


03 


46.0 


0.1659 


± 





0079 


4 


OE- 


03 


46.0 


0.1339 


± 





0084 


5 


9E- 


03 


46.0 


0.1127 


± 





0086 


1 


OE- 


02 


46.0 


0.1026 


± 





0058 


5 


9E 


03 


130.0 


0.140 


± 





015 


1 


OE- 


02 


130.0 


0.1240 


± 





0086 


2 


5E- 


02 


130.0 


0.0880 


± 





0071 
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Table 12: The structure function F 2 as a function of Xl in bins of x and 
Q 2 (DIS sample), for protons in a restricted p\ range, p\ < 0.04 GeV 2 . The two 
rightmost values indicate the statistical and systematic uncertainties, respectively. 



Q 2 (GeV 2 ) 


X 


XL 








F 


T P f 3 ) 
2 








4.4 


3 


3E 


04 





73 





0675 


± 





0080 


± 





0157 


4.4 


3 


3E 


04 





78 





0724 


± 





0078 


± 





0116 


4.4 


3 


3E 


04 





83 





0756 


± 





0082 


± 





0106 


4.4 


3 


3E 


04 





88 





0585 


± 





0082 


± 





0098 


4.4 


1 


0E 


03 





73 





0489 


± 





0077 


± 





0125 


4.4 


1 


OE 


03 





78 





0552 


± 





0077 


± 





0101 


4.4 


1 


OE 


03 





83 





0605 


± 





0083 


± 





0099 


4.4 


1 


OE 


03 





88 





0460 


± 





0083 


± 





0093 


7.5 


3 


3E 


04 





73 





0962 


± 





0079 


± 





0209 


7.5 


3 


3E 


04 





78 





0955 


± 





0074 


± 





0136 


7.5 


3 


3E 


04 





83 





0822 


± 





0070 


± 





0101 


7.5 


3 


3E 


04 





88 





0855 


± 





0083 


± 





0115 


7.5 


1 


OE 


03 





73 





0899 


± 





0085 


± 





0200 


7.5 


1 


OE 


03 





78 





0737 


± 





0072 


± 





0114 


7.5 


1 


OE 


03 





83 





0682 


± 





0071 


± 





0094 


7.5 


1 


OE 


03 





88 





0715 


± 





0084 


± 





0107 


7.5 


3 


3E 


03 





73 





0587 


± 





0081 


± 





0143 


7.5 


3 


3E 


03 





78 





0563 


± 





0075 


± 





0100 


7.5 


3 


3E 


03 





83 





0505 


± 





0072 


± 





0085 


7.5 


3 


3E 


03 





88 





0619 


± 





0093 


± 





0109 


13.3 


1 


OE 


03 





73 





0957 


± 





0078 


± 





0208 


13.3 


1 


OE 


03 





78 





1023 


± 





0076 


± 





0144 


13.3 


1 


OE 


03 





83 





0886 


± 





0072 


± 





0107 


13.3 


1 


OE 


03 





88 





0729 


± 





0074 


± 





0101 


13.3 


3 


3E 


03 





73 





0666 


± 





0061 


± 





0147 


13.3 


3 


3E 


03 





78 





0698 


± 





0059 


± 





0102 


13.3 


3 


3E 


03 





83 





0677 


± 





0059 


± 





0084 


13.3 


3 


3E 


03 





88 





0658 


± 





0067 


± 





0091 


28.6 


1 


OE 


03 





73 





1249 


± 





0125 


± 





0281 


28.6 


1 


OE 


03 





78 





1290 


± 





0119 


± 





0195 


28.6 


1 


OE 


03 





83 





1007 


± 





0107 


± 





0140 


28.6 


1 


OE 


03 





88 





0956 


± 





0121 


± 





0150 


28.6 


3 


3E 


03 





73 





0906 


± 





0089 


± 





0203 


28.6 


3 


3E 


03 





78 





0934 


± 





0085 


± 





0140 


28.6 


3 


3E 


03 





83 





0701 


± 





0075 


± 





0097 


28.6 


3 


3E 


03 





88 





0758 


± 





0089 


± 





0114 
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Table 13: Fraction of leading-proton DIS events with exactly two jets with E T > 
4 GeV, rfp, as a function of xl for p\ < 0.5 GeV 2 . The two rightmost values 
indicate the statistical and systematic uncertainties, respectively. The systematic 
uncertainties are highly correlated. 



X L 


iet 
r LP 


0.645 


0.0209 ± 0.0031 ± 0.0053 


0.735 


0.0254 ± 0.0023 ± 0.0062 


0.825 


0.0242 ± 0.0019 ± 0.0060 


0.920 


0.0196 ± 0.0026 ± 0.0050 



Table 14: Fraction of leading-proton DIS events with exactly two jets with Et > 
4 GeV, 7" j^p j (IS CL function of p\ for 0.6 < xl < 0.97. The two rightmost values 
indicate the statistical and systematic uncertainties, respectively. The systematic 
uncertainties are highly correlated. 



p\ (GeV 2 ) 


jet 
r LP 


0.0105 
0.0355 
0.0900 
0.3150 


0.0201 ± 0.0020 ± 0.0037 
0.0204 ± 0.0022 ± 0.0064 
0.0259 ± 0.0024 ± 0.0057 
0.0350 ± 0.0032 ± 0.0083 
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Table 15: Ratio of the yield of DIS events with exactly two jets with E T > 4 GeV 
and an LPS proton to the yield of DIS events with exactly two jets, also with 
Et > 4 GeV, rj^, as a function of Et of the higher-energy jet. The statistical 
uncertainty is given. A fully correlated systematic uncertainty of ±13% is not 
included. 



E T (GeV) 


r LP 
'jet 


4.8 


0.126 ± 0.018 


5.8 


0.122 ± 0.017 


6.8 


0.130 ± 0.019 


8.0 


0.124 ± 0.018 


12.3 


0.088 ± 0.014 



Table 16: Ratio of the yield of DIS events with exactly two jets with Et > 4 GeV 
and an LPS proton to the yield of DIS events with exactly two jets, also with 
E T > 4 GeV, rj^, as a function of Q 2 . The statistical uncertainty is given. A fully 
correlated systematic uncertainty of ±13% is not included. 



Q 2 (GeV 2 ) 


/>! LP 
'jet 


6.6 


0.123 ± 0.013 


19.4 


0.097 ± 0.017 


36.0 


0.111 ± 0.017 


106.8 


0.121 ± 0.015 



Table 17: Ratio of the yield of DIS events with exactly two jets with Et > 4 GeV 
and an LPS proton to the yield of DIS events with exactly two jets, also with 
Et > 4 GeV, rj^, as a function of x. The statistical uncertainty is given. A fully 
correlated systematic uncertainty of ±13% is not included. 



X 


,ViLP 

'jet 


0.00027 


0.131 ± 0.014 


0.00093 


0.109 ± 0.015 


0.0022 


0.107 ± 0.015 


0.0079 


0.109 ± 0.017 
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Figure 1: Schematic diagram of the reaction e + p — > e + Xp. 
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Figure 2: Distributions of the variables a) E' e , b) $ e and c) y,jB for the recon- 
structed BPC data (squares) and the simulated events (EPSOFT), shown as the 
shaded histograms (normalised to the data); d) E' e , e) ~d e and f) 7^ for the re- 
constructed DIS data (dots) and the simulated events (RAPGAP), shown as the 
hatched histograms (normalised to the data). 
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Figure 3: Distributions of the variables a) xl, b) p%, c) Ap ipc and d) A p i ane , 
for the reconstructed BPC data (squares) and for the simulated events (EPSOFT) 
shown as the shaded histogram (normalised to the data). The arrows indicate the 
minimum allowed values of A pipe and A p i ane (see text). 
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Figure 4: LPS geometrical acceptance for different Xl values as a function of 
Px and Py . The shaded areas indicate the regions of acceptance. The dashed 
circles indicate the limits of the Pt bins used in the analysis (the bin edges are 
0.124 GeV, 0.199 GeV, 0.280 GeV, 0.370 GeV, 0.470 GeV, 0.507 GeV, 0.707 GeV). 
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Figure 5: a) Distribution of the DIS events in the ^ max , xl) plane. The horizontal 
and vertical lines indicate ?7 max = 2.5 and xl = 0.97 ; respectively, b) Fraction of 
events with r] max < 2.5 for both the BPC and the DIS samples. The BPC data are 
slightly shifted in Xl for clarity of presentation. The bars indicate the statistical 
uncertainties; systematic uncertainties mostly cancel in the ratio. 
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Figure 6: The normalised cross-section (l/er tot ) • do^* p _,x-pl 'dxi for the BPC and 
DIS data compared to the pp data [54] in the region p\ < 0.5 GeV 2 . The inner 
bars indicate the statistical uncertainties and the outer bars are the statistical and 
systematic uncertainties summed in quadrature. 
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Figure 7: The normalised cross-section (l/cr tot ) • da^* p ^xp/dxL for the pho- 
toproduction, BPC and DIS data compared to the pp data [54] in the region 
p\ < 0.04 GeV 2 . The inner bars indicate the statistical uncertainties and the 
outer bars are the statistical and systematic uncertainties summed in quadrature. 
The HI results [7] are also shown. 
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Figure 8: The normalised cross-section (l/cr tot ) • da 1 * p ^ i xp/d x L for the DIS data 
(as shown in Fig. 6) compared to a) the model of Duraes et al. (solid curve), 
LEPT06.5 (dashed curve) and DJANGO (dot-dashed curve), and b) to the model 
of Szczurek et al. (solid curve). For the latter, the individual contributions of 
Pomeron, Reggeon and pion exchanges are indicated; for pion exchange, the con- 
tribution of final states with isospin 1=1/2 and 1 = 3/2 are shown separately. 
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Figure 9: The differential cross-section da^p^xp/dPr f or several xl bins for 
the BPC sample. The lines represent the results of fits to the functional form 
do-y*p^xp/ dPr e~ bpT . The fitted values of b and their statistical uncertainties 
are also given. The inner bars indicate the size of the statistical uncertainties, the 
outer bars show the statistical and systematic uncertainties summed in quadrature. 
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Figure 10: The differential cross-section day*p->xp/ 1 dVr for several xl bins for 
the DIS sample. The lines represent the results of fits to the functional form 
da^p^xp/ dPr ^ e~ bpT . The fitted values of b and their statistical uncertainties 
are also given. The inner bars indicate the size of the statistical uncertainties, the 
outer bars show the statistical and systematic uncertainties summed in quadrature. 
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Figure 12: Tfte ratao r LP ( 3 ) = / F 2 as a function of Xl in bins of x and Q 2 

(BPC sample), for protons withp\ < 0.5 GeV 2 . The inner bars show the statistical 
uncertainties and the outer bars the statistical and systematic uncertainties added 
in quadrature. The dashed line r LP ^ = 0.4 is overlaid to guide the eye. 
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Figure 14: The ratio r LP ^ = / F 2 as a function of x for fixed Q 2 values, 

for protons with 0.6 < xl < 0.97 and p?r < 0.5 GeV 2 . The error bars show the 
statistical uncertainties. A fully correlated systematic uncertainty of ±13% is not 
shown. The horizontal lines r LF ^ = 0.10 and r LF ^ = 0.15 are overlaid to guide 
the eye. 
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Figure 15: The average ratio (r LP ( 2 )) = F 2 LP( ^ / F 2 as a function ofQ 2 . The error 
bars show the statistical uncertainties. A fully correlated systematic uncertainty of 
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Figure 16: The structure-function F 2 as a function of x for 0.6 < xl < 0.97 
and p\ < 0.5 GeV 2 . The bands show the one-standard-deviation limits of the 
F 2 parametrisations used, scaled by the average value of r LF ^ ^r LF ^) ~ 0.13). 
The error bars show the statistical uncertainties. A fully correlated systematic 
uncertainty of ±13% is not shown. 
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Figure 17: The structure function F 2 LP( ' 3 ^ as a function of Xl in bins of x and 
Q 2 (DIS sample), for protons in a restricted p\ range, p\ < 0.04 GeV 2 . The inner 
bars show the statistical uncertainties and the outer bars are the statistical and 
systematic uncertainties added in quadrature. The HI results [5] are also shown. 
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Figure 18: Fraction of leading-proton DIS events with exactly two jets with 
E T > 4 GeV, 7lp, as a function of a) xl for -p\ < 0.5 GeV 2 and b) for 
0.6 < xl < 0.97. The inner bars show the statistical uncertainties and the outer 
bars show the statistical and systematic uncertainties added in quadrature. The 
systematic uncertainties are highly correlated. 
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Figure 19: Ratio of the yield of DIS events with exactly two jets with Et > 4 GeV 
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Ex > 4 GeV, rj^, as a function of a) E? of the higher- energy jet, b) Q 2 and c) 
x. The error bars show the statistical uncertainties. A fully correlated systematic 
uncertainty of ±13% is not shown. 
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